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Abstract 

Path-integral Monte Carlo calculations have been performed to study the 4 He adsorption on 
both sides of a single a-graphyne sheet. For investigation of the interlayer correlation between 
the upper and the lower monolayer of 4 He adatoms, the 4 He-substrate interaction is described 
by the sum of the 4 He-C interatomic pair potentials, for which we use both Lennard-Jones and 
Yukawa-6 anisotropic potentials. When the lower 4 He layer is a C 4 / 3 commensurate solid, the 
upper-layer 4 He atoms are found to form a Kagome lattice structure at a Mott insulating density 
of 0.0706 A' 2 , and a commensurate solid at an areal density of 0.0941 A -2 for both substrate 
potentials. The correlation between upper- and lower-layer pseudospins, which were introduced 
in Ref. [1] for two degenerate configurations of three 4 He atoms in a hexagonal cell, depends on 
the substrate potential used; With the substrate potential based on the anisotropic Yukawa-6 pair 
potentials, the Ising pseudo-spins of both 4 He layers are found to be anti-parallel to each other 
while the parallel and anti-parallel pseudo-spin alignments between the two 4 He layers are nearly 
degenerate with the Lennard-Jones potentials. This is attributed to the difference in the interlayer 
distance, which is ~ 4 A with the Yukawa-6 substrate potential but as large as ~ 4.8 A with the 
Lennard-Jones potential. 
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I. INTRODUCTION 


Among many substrates, graphite has long served as a test bed to investigate low¬ 
dimensional quantum fluids because of its strong binding of adsorbates. Up to seven 


distinct 4 He layers were observed on grap 
a quasi-two-dimensional quantum system 


ihite and each helium layer is considered to be 
2J]. The first 4 He adlayer on graphite shows a 
commensurate-incommensurate solid transition as the helium coverage increases Bfl. Re¬ 
cently, a series of theoretical calculations have been performed to study the 4 He adsorption 
on newly-synthesized (or -proposed) low-dimensional carbon substrates such as graphene jfj- 
8 ], graphynes jl, 9], carbon nanotubes jlO, llj], and fullerene molecules 12 15]. The phase 
diagrams of the 4 He layers adsorbed on graphene were predicted to be very similar to those 
of the corresponding layers on graphite; the monolayer of 4 He adatoms shows a C 1/3 com¬ 
mensurate structure at the areal density of 0.0636 A ~ 2 and goes through various domain-wall 
phases before crystallizing into an incommensnrate triangular solid near its completion flfl. 

Graphyne is a two-dimensional (2D) network of sp- and sp 2 -bonded C atoms [16l . Il7 ] 
which could be permeable to a 4 He gas unlike graphene. Despite much experimental effort 
motivated by some promising theoretical predictions for graphyne as new Dirac materi¬ 


als 


181-3 and high-capacity energy storage materials (2 ll l23|. there has been no successful 


report yet for fabrication of extended 2D graphynes. However, some flakes or building blocks 
of finite-size graphynes have been synthesized 24-26), leading to a belief that graphynes will 
be fabricated in the near future. On the surface of 7 -graphyne, which is the most stable 
graphyne structure according to quantum Monte Carlo calculations 27), the 4 He monolayer 
was predicted to exhibit a richer phase diagram than the corresponding layer on graphene or 
graphite, including various commensurate and incommensurate structures depending on the 
helium density 9]. Recently one of us performed path-integral Monte Carlo (PIMC) calcula¬ 


tions for 4 He atoms adsorbed on a AB-stacked bilayer a-graphyne jl| , which is a hybridized 
honeycomb structure with each hexagon side consisting of one sp 2 and two sp C atoms. 
It was found that the 4 He monolayer was in a Mott insulating state at an areal density of 
0.0706 A -2 while a commensurate solid was realized at 0.0941 A~ 2 . Introducing Ising pseudo¬ 
spin degrees of freedom for two degenerate configurations for three 4 He atoms occupying a 
hexagonal cell (see Fig. 3 of Ref. |l]), this Mott-insulator to commensurate-solid transition 
was interpreted as a symmetry breaking process from a spin liquid of geometrically-frustrated 
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Q. 


antiferromagnets to a spin-aligned ferromagnet 

One interesting feature of a 2D carbon structure, such as graphene and graphyne, is that 

and 4 ffe atoms can be coated on both sides. Noting that 


it can be suspended in the air 


some new physics could emerge as a result of interlayer correlation between opposite-side 
4 ffe layers, some theoretical studies were recently done for the 4 He adsorption on both sides 
of a single graphene sheet. Markic et al. found that the correlation between two 4 He clusters 
adsorbed on opposite sides of graphene, 5 ~ 6 A apart from each other, was quite weak 
as evidenced by peakless pair distribution functions 29]. A weak correlation between two 
4 ffe systems on the opposite sides of graphene was also predicted by Gordillo’s diffusion 
Monte Carlo calculations, which showed that the phase diagram of the 4 He monolayer on 
graphene would not be affected by the 4 He adsorption on the other side |^]. In this paper 
we report PIMC study of the 4 He adsorption on both sides of a single a-graphyne sheet. 
Because a-graphyne is more porous than graphene, 4 He atoms can penetrate through gra¬ 
phyne to allow physical exchanges among 4 He atoms on opposite sides. This could result in 
stronger interlayer correlation than the corresponding systems on graphene. We find that 
4 He atoms in a Mott-insulating state form a 2D Kagome lattice as a result of the interlayer 
correlation when the opposite-side 4 He layer is a C 4/3 commensurate solid, a ferromagnetic 
state in a pseudospin terminology. Effects of the interlayer correlation between two ferro¬ 
magnetic C 4/3 solids are found to depend on the substrate potential used; the parallel and 
the antiparallel pseudospin alignments between two 4 He layers are nearly degenerate with 
the substrate potential based on the Lennard-Jones (LJ) 4 He-C pair potentials while the 
antiparallel alignment is favored with the one described by the Yukawa-6 pair potentials. 
The vacancy formation in a 4 He layer on a-graphyne is also found to be affected by the 
presence of the opposite-side 4 He layer. 

In the following section, we outline our approach and some computational details. The 
PIMC results along with the related discussions are presented in detail in Sec. IIII1 We 
summarize our Endings in Sec. IIV1 


II. METHODOLOGY 

In this study, a single a-graphyne sheet is set to be at z = 0. The 4 He-graphyne interaction 
is assumed to be a pairwise sum of interatomic potentials between the carbon atoms and 
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a 4 He atom, which has been widely used to describe the interaction between a 4 He atom 


and a carbon substrate 


, 7, 9, 29, 3o|. For the 4 He- 


C interatomic pair potential, we 


employ two anisotropic potentials proposed by Carlos and Cole 31, [32j, i.e., a 6-12 LJ 


potential and a Yukawa-6 potential. For the computational convenience our previous study 
for the 4 He monolayer on bilayer a-graphyne was done with only isotropic parts of the LJ 
pair potential. However, the original interatomic pair potentials of Carlos and Cole include 
anisotropic parts to fit helium scattering data from graphite surfaces. Even though the 
inclusion of the anisotropic parts of the interatomic potentials has little effect on quantum 
phases displayed by the 4 He layer on one side of a-graphyne, it allows some 4 He atoms to 
be closer to the substrate, resulting in stronger correlation between two 4 He layers on the 
opposite sides (the minima of the substrate potential made of the anisotropic pair potentials 
are deeper and closer to graphyne than the corresponding ones based on only isotropic parts 
of the pair potentials). This leads to our decision of using the substrate potentials based 
on fully-anisotropic interatomic pair potentials, which should give a better description of 
the interlayer correlation. Furthermore, since the LJ and the Yukawa-6 potentials used in 
this study were based on an interaction between helium and sp 2 -bonded carbon atoms in 
graphite, we tested the sensitivity of our modelling of 4 He-graphyne potentials to the well 
depth of the pair potentials. Although decrease in the well depth yields more fluctuations 
in 4 He density distributions, the density modulations are found to change only little and 
our main results presented below are still, at least qualitatively, valid. For the 4 He- 4 He 
interaction, we use a well-known Aziz potential 33]. 

In the discrete path-integral representation, the thermal density matrix at a low tem¬ 
perature is expressed by a convolution of M high-temperature density matrices with an 
imaginary time step of r = 1/(MksT) 34]. While the isotropic parts of 4 He-C pair po¬ 
tentials along with the 4 He- 4 He potential pair potentials are used to compute the exact 


two-body density matrices |34l . [35] at the high temperature MT, their anisotropic parts are 
treated with the primitive approximation 34]]. This is found to give accurate description 
of both 4 He- 4 He and 4 He-graphyne interactions with a time step of (tA;#) -1 = 80 K. We 
employ the multilevel Metropolis algorithm to sample the imaginary time paths along with 
permutations among 4 He atoms as described in Ref. 34]. To minimize finite size effects, 
periodic boundary conditions are applied along the lateral directions. 
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III. RESULTS 


The PIMC calculations were done with a fixed 3x2 simulation cell with dimensions of 
21.01 x 24.26 A 2 , the same as in our previous study for 4 He on bilayer a-graphyne [ 1 ]. We 
focus on the interlayer correlation between two 4 He layers on the opposite sides of graphyne, 
which are either in a Mott-insulating state or a pseudospin-aligned commensurate solid 
state. The results obtained with two different substrate potentials, the LJ potential and the 
Yukawa -6 one, are presented separately below. 

A. Lennard-Jones substrate potential 

For PIMC calculations with the LJ substrate potential, we first prepare the a-graphyne 
surface whose in-plane hexagon center is occupied by a single 4 He atom and whose lower 
side is coated with a monolayer of 4 He atoms constituting a C 4/3 commensurate solid while 
each of the in-plane centers is occupied by a single 4 He atom. The simulations for the 4 He 
adsorptions on the upper side of the prepared graphyne surface begin from an initial config¬ 
uration of 4 He atoms being randomly distributed at the distances far away from graphyne. 
Figure |Tj presents one-dimensional (ID) density distributions of 4 He atoms, as a function 
of the vertical coordinate z along the direction perpendicular to the graphyne surface, for 
two different combinations of particle numbers per simulation cell. Two distinct density 
peaks, which correspond to the upper and the lower 4 He layers, are observed on the oppo¬ 
site sides of graphyne. Note that 36 and 48 4 He atoms per simulation cell correspond to the 
Mott-insulating density of 0.0706 A ~ 2 and the C 4 / 3 commensurate density of 0.0941 A” 2 , 
respectively. The additional density peak at z = 0 corresponds to the zeroth layer consisting 
of 4 He atoms embedded onto the in-plane hexagon centers, which was also observed in our 
previous study for 4 He adatoms on a bilayer a-graphyne g Since the p^o-pea. distance 
between the upper and the lower layers is estimated to be about 4.8 A, one can expect that 
the van der Waals interaction between 4 He atoms on the opposite sides is weakly attractive 
(note that the Aziz potential we used for the 4 He- 4 He interaction has a minimum value at 
r ~ 3.0 A). In addition, the clear separation between the adjacent density peaks in Figure [Tj 
suggests that exchange couplings among 4 He atoms in different layers are nearly absent and 
any correlation between the upper and the lower layers, if it exists, should stem mostly from 
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the weakly-attractive 4 He- 4 He interaction rather than particle exchanges. 

Figure [2] shows two-dimensional (2D) density distributions of the upper-layer 4 He atoms, 
while the lower-layer densrty peaks represented by the white stars constitute a C 4/3 com- 
mensurate structure with all pseudospins being in the spin-up state (see Fig. 3 of Ref. [Ij). 
ffere a distinct density peak in each plot represents an occupancy of a single 4 He atom. 
At an areal density of 0.0706 A -2 , every hexagonal cell of graphyne is seen in Fig. [21(a) to 
accommodate three upper-layer 4 He atoms, which is a manifestation of a Mott-insulating 
state. Without the lower 4 He layer in a pseudospin-aligned commensurate solid state, this 
Mott-insulating state is a nonmagnetic spin liquid of frustrated antiferromagnets in terms 
of pseudospin degrees of freedom (see Fig. 2(d) in Ref. [1]). ffowever, in the presence of 
the pseudospin-aligned lower 4 He layer, the upper-layer pseudospins are shown in Fig. [21(a) 
to be aligned in the same direction as the lower-layer ones. Our PIMC simulations at 
T = 0.5 K have also produced the antiparallel pseudospin alignment between the two 4 He 
layers. This is understood by the fact that the parallel alignment is energetically favored 
only by ~ 0.3 K per an upper-layer helium atom over the antiparallel alignment, i.e., two 
pseudospin alignments are nearly degenerate. We here note that the upper-layer 4 He atoms 
in a pseudospin-aligned Mott insulating state of Fig. [21(a) constitute a 2D Kagome lattice. 
This is also true when the upper-layer pseudospins are aligned in the opposite direction to 
the lower-layer ones. Therefore one can conclude that as a result of the interlayer correlation, 
the upper-layer 4 He atoms form a Kagome lattice structure at the Mott insulating density of 
0.0706 A -2 when the lower 4 He atoms constitute a pseudospin-aligned C 4/3 commensurate 
solid. 

The interlayer correlation between two ferromagnetic C 4/3 commensurate solids is also 
analyzed. Figure [21(b) presents the 2D density distribution of the upper-layer 4 He atoms 
at the areal density of 0.0941 A^ 2 , where they constitute a 4/3 commensurate solid. The 
upper-layer pseudospins are seen to be aligned in the same direction as the lower-layer ones. 
Similarly to the case of the Mott-insulating state, the antiparallel pseudospin alignment was 
also observed in our simulations. These two pseudospin alignments are more degenerate 
(the parallel alignment was found to be preferred by ~ 0.11 K per upper-layer 4 He atom), 
than the parallel and antiparallcl pseudospin alignments between a Mott-insulator and a 
C 4/3 solid in Fig. [21(a). This can be understood by the fact that unlike 4 He atoms inside a 
hexagonal cell, upper-layer 4 He atoms at the vertices of the graphyne hexagons in a C 4/3 
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solid state prefer the other sublattice sites over the ones occupied by the corresponding 
lower-layer helium atoms. We note that once a pseudospin alignment between the two 4/3 
commensurate 4 He solids is established, either parallel or antiparallel to each other, energy 
barrier is too large to reverse one alignment to the other. 

Now we analyze the effects of the interlayer correlation on the formation of vacancies 
especially in a Mott insulator with the Kagome lattice structure. Figure [3] shows the 2D 
density distribution of 35 upper-layer 4 He atoms, one less than the Mott-insulating case, per 
3x2 simulation cell. So the upper-layer Mott insulator contains one vacancy per simulation 
cell while the lower 4 He layer is the same 4/3 commensurate solid as in Fig. [2fa). One 
can see that every hexagonal cell, except one, is seen to accommodate three upper-layer 4 He 
atoms and its pseudospin is aligned in the same direction as those of the ferromagnetic lower 
layer. As shown in Fig. |3j one cell involving only two upper-layer atoms does not show the 
clear pseudospin alignment. This tells us that the Kagome lattice structure is sustained even 
with the creation of vacancies but those isolated vacancies are restricted at one triangle of 
this trihexagonal tiling structure without hopping to the neighboring sites because of high 
potential barrier provided by the graphyne surface. 


B. Yukawa-6 substrate potential 


Our PIMC simulations with the Yukawa-6 substrate potential start from an initial con¬ 
figuration of N up and N dn 4 He atoms being distributed randomly on the upper and the lower 
side of a-graphyne. respectively. Figure [4] presents the ID 1 He density distributions as a 
function of the vertical coordinate z for two different values of N up while N dn is fixed to 
48 per 3x2 simulation cell. Unlike Fig. [1] for the LJ substrate potential, only two density 
peaks on the opposite sides of graphyne are observed without the zeroth layer consisting 
of 4 He atoms embedded onto the in-plane hexagon centers. The absence of the zeroth 4 He 
layer is attributed to the fact that the Yukawa-6 substrate potential is more slowly varying 
near the potential minima, that is, the hexagon centers, than the LJ substrate potential. 
Note that the Yukawa-6 4 He-C pair potential is less repulsive at short distances than the 
LJ interatomic pair potential 32], The distance between two density peaks is ~ 4.0 A, 
shorter than the corresponding distance for the LJ substrate potential. Furthermore, there 
is significant overlap between the two density peaks that are broader than those in Fig. Q] 
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This indicates large quantum fluctuations of 4 He adatoms along the vertical direction, which 
could result in frequent particle exchanges between these two layers. 

Figure 0 shows 2D density distributions of the upper-layer 4 He atoms, while the lower- 
layer density peaks are represented by the white stars. Even with the Yukawa -6 substrate 
potential, the lower 4 He layer is seen to constitute a pseudospin-aligned a C 4 / 3 commensurate 
solid at the areal density of 0.0941 A~ 2 . This provides another confirmation to the conclusion 
of Ref. jl] that most of quantum phases manifested in the 4 He monolayer on the Q-graphyne 
surface, such as a Mott insulator, commensurate solids, and pseudospin degrees of freedom, 
are not sensitive to the specifics of the substrate potential but are determined mostly by the 
surface geometry. It is also shown in Fig. |5Ja) that the upper-layer 4 He atoms under the 
Yukawa -6 substrate potential are in a Mott-insulating state at the areal density of 0.0706 A -2 
with each hexagonal cell accommodating three 4 He atoms. Unlike Fig. [2fa), however, all 
upper-layer pseudospins are aligned in the opposite direction to those of the ferromagnetic 
lower-layer commensurate solid. We understand that the increase in the effective hard-core 
radii of 4 He adatoms due to larger quantum fluctuations, along with a shorter interlayer 
distance, causes the antiparallel pseudospin alignment to be favored under the Yukawa- 
6 substrate potential. As observed with the LJ substrate potential, the upper-layer 4 He 
adatoms in the pseudospin-aligned Mott insulating state constitute a 2D Kagome lattice 
structure. The interlayer correlation that favors the antiparallel pseudospin alignment is 
more evident between two ferromagnetic C 4/3 solids. We observe in Fig. [S^b) that the 
pseudospins of a upper-layer 4/3 commensurate solid are in spin-down state while the lower- 
layer pseudospins are in spin-up state. This antiparallel pseudospin alignment between the 
two 4 He adlayers corresponds to the AB stacking of two triangular solids. So we conclude 
that with the Yukawa -6 substrate potential, the AB stacking is preferred to the AA stacking 
between two C 4/3 triangular 4 He solids while these two stacking orders are nearly degenerate 
with the LJ potential as discussed in Sec. IIIIA1 

We now try to create a single vacancy in the upper-layer Mott insulator, i.e., the Kagome 
lattice structure, by putting only 35 4 He atoms on the upper side in our simulation. Figure [ 6 ] 
presents the 2D density distribution of the upper-layer 4 He atoms along with the lower-layer 
density peaks of a pseudospin-up C 4/3 commensurate solid. It is shown that there are still 
three upper-layer density peaks per every hexagonal cell, corresponding to a configuration 
of a Mott insulator with all pseudospins being in the down state. However, one lower-layer 





density peak or one white star is missing at a vertex of the graphyne honeycomb structure 
(see the bottom right corner), which prevents the lower-layer 4 He atoms from forming a 
perfect C 4 / 3 triangular lattice. This suggests that when a single vacancy is created in a 
upper-layer Kagome lattice, one lower-layer 'lie atom moves to the upper layer to form a 
perfect upper-layer lattice structure while a localized vacancy is created in the lower-layer 
4/3 commensurate solid. This lower-layer vacancy is found at a vertex site on top of a carbon 
atom because it is a less favorable site for a 4 He adatom than a site inside a hexagonal cell. 
The layer-to-layer hopping of a vacancy can be understood by a short interlayer distance 
and large quantum fluctuations along the vertical direction under the Yukawa -6 substrate 
potential. 


IV. CONCLUSION 


According to our PIMC calculations of using two different 4 He-substrate potentials, 4 He 
atoms form distinct layers on both sides of a single o-graphyne sheet. Regardless of the 
substrate potential used, the upper-layer 4 He atoms form a 2D Kagome structure at the 
Mott-insulating density of 0.0706 A -2 as a result of the interlayer correlation when the lower 
layer is a pseudospin-aligned C 4 / 3 commensurate solid. Since the interaction of 3 He atoms 
with a substrate or between themselves is similar to the corresponding interaction for 4 He, 
the same Kagome lattice structure is expected to be formed in the fermionic counterpart of 
the upper 4 He layer , i.e., a 4 He upper layer adsorbed on a-graphyne, when its lower side is 
coated with the C 4 / 3 commensurate helium solid. We speculate that some novel phenomena 
related with a geometrically-frustrated antiferromagnetism such as quantum spin liquids 36], 
could emerge in this 3 He Kagome lattice. 

The interlayer correlation results in different stacking orders between two C 4 / 3 commensu¬ 
rate triangular solids on the opposite sides of graphyne, depending on the substrate potential; 
with the Yukawa -6 potential, the AB stacking (an antiparallcl pseudospin alignment between 
two 4 He solids) is found to be favored but both AA (a parallel pseudospin alignment) and 
AB stacking configurations are nearly degenerate with the LJ substrate potential. This is 
attributed to the difference between two substrate potentials in the interlayer distance as 
well as in the magnitude of quantum fluctuations along the vertical direction. A more ac¬ 
curate 4 He-graphyne potential would be required to draw a definite conclusion about the 
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preferred stacking order of two commensurate triangular 4 He solids on the opposite sides of 
a-graphyne. 

Recent theoretical studies done by Markic et al. 29] and Gordillo 30] reported that the 
interlayer correlation between two 4 He systems adsorbed on both sides of graphene was very 
weak and quantum phase diagram of one 4 He layer would not be affected by the presence 
of the opposite-side 4 He layer. On the other hand, our PIMC calculations have revealed 
some significant effects of the interlayer correlation on the structural properties of the 4 He 
monolayers on a-graphyne. This difference is understood to be due to much more porous 
nature of a-graphyne than graphene. 


Acknowledgments 

This work was supported by Konkuk University. We also acknowledge the support from 
the Supercomputing Center/Korea Institute of Science and Technology Information with 
supercomputing resources including technical support (KSC-2013-C3-033). 


[1] Y. Kwon, H. Shin, and H. Lee, Phys. Rev. B 88, 201403(R) (2013). 

[2] G. Zimmerli, G. Mistura, and M. H. W. Chan, Phys. Rev. Lett. 68, 60 (1992). 

[3] D. S. Greywall and P. A. Busch, Phys. Rev. Lett. 67, 3535 (1991). 

[4] D. S. Greywall, Phys. Rev. B 47, 309 (1993). 

[5] P. A. Crowell and J. D. Reppy, Phys. Rev. B 53, 2701 (1996). 

[6] M. C. Gordillo and J. Boronat, Phys. Rev. Lett. 102, 085303 (2009). 

[7] Y. Kwon and D. M. Ceperley, Phys. Rev. B 85, 224501 (2012). 

[8] J. Happacher, P. Corboz, M. Boninsegni, and L. Pollet, Phys. Rev. B 87, 094514 (2013). 

[9] J. Ahn, H. Lee, and Y. Kwon, Phys. Rev. B 90, 075433 (2014). 

[10] M. W. Cole, V. H. Crespi, G. Stan, C. Ebner, J. M. Hartman, S. Moroni, and M. Boninsegni, 
Phys. Rev. Lett. 84, 3883 (2000). 

[11] M. C. Gordillo, Phys. Rev. Lett. 101, 046102 (2008). 

[12] Y. Kwon and H. Shin, Phys. Rev. B 82, 172506 (2010). 

[13] H. Shin and Y. Kwon, J. Chem. Phys. 136, 064514 (2012). 


10 





[14] B. Kim and Y. Kwon, J. Low Temp. Phys. 171, 599 (2013). 

[15] S. Park and Y. Kwon, Phys. Rev. E. 89, 042118 (2014). 

[16] R. H. Baughman, H. Eckhardt, and M. Kertesz, J. Chem. Phys. 87, 6687 (1987). 

[17] V. R. Coluci, S. F. Braga, S. B. Legoas, D. S. Galvao, and R. H. Baughman, Nanotechnology 
15, S142 (2004). 

[18] D. Malko, C. Neiss, F. Vines, and A. Gorling, Phys. Rev. Lett. 108, 086804 (2012). 

[19] B. G. Kim and H. J. Choi, Phys. Rev. B 86, 115435 (2012). 

[20] J. Chen, J. Xi, D. Wang, and Z. Shuai, J. Phys. Chem. Lett. 4, 1443 (2013). 

[21] H. Zhang, M. Zhao, X. He, Z. Wang, X. Zhang, and X. Liu, J. Phys. Chem. C 115, 8845 

( 2011 ). 

[22] H. Hwang, Y. Kwon, and H. Lee, J. Phys. Chem. C 116, 20220 (2012). 

[23] H. Hwang, J. Koo, M. Park, N. Park, Y. Kwon, and H. Lee, J. Phys. Chem. C 117, 6919 
(2013). 

[24] Y. Tobe, I. Ohki, M. Sonoda, H. Niino, T. Sato, and T. Wakabayashi, J. Am. Chem. Soc. 
125, 5614 (2003). 

[25] M. M. Haley, Pure Appl. Chem. 80, 519 (2008). 

[26] F. Diederich and M. Kivala, Adv. Mater. 22, 803 (2010). 

[27] H. Shin, S. Kang, J. Koo, H. Lee, J. Kim, and Y. Kwon, J. Chem. Phys. 140, 114702 (2014). 

[28] J. C. Meyer, A. K. Geirn, M. I. Katsnelson, K. S. Novoselov, T. J. Booth, and S. Roth, Nature 
446, 60 (2007). 

[29] L. Vranjes Markic, P. Stipanovic, I. Beslic, and R. E. Zillich, Phys. Rev. B 88, 125416 (2013). 

[30] M. C. Gordillo, Phys. Rev. B 89, 155401 (2014). 

[31] W. E. Carlos and M. W. Cole, Phys. Rev. Lett. 43, 697 (1979). 

[32] W. E. Carlos and M. W. Cole, Surf. Sci. 91, 339 (1980). 

[33] R. A. Aziz, M. J. Slarnan, A. Koide, A. R. Allnatt, and W. J. Meath, Mol. Phys. 77, 321 
(1992). 

[34] D. M. Ceperley, Rev. Mod. Phys. 67, 279 (1995). 

[35] R. E. Zillich, F. Paesani, Y. Kwon, and K. B. Whaley, J. Chem. Phys. 123, 114301 (2005). 

[36] L. Balents, Nature 464, 199 (2010). 


11 




Z[A] 


FIG. 1: (Color online) One-dimensional density distributions of 4 He atoms adsorbed on both sides 
of a-graphyne as a function of the vertical coordinate z perpendicular to the graphyne surface, 
which were computed with the LJ substrate potential. Here IV U p and N^ n represent the number of 
4 He atoms per 3x2 simulation cell in the upper and the lower 4 He layer, respectively. Additional 
12 4 He atoms per simulation cell are involved to form the zeroth layer around z = 0 where one 4 He 
atom is embedded at every hexagon center (see Ref. [11). 
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FIG. 2: (Color online) Contour plots of two-dimensional density distributions of 4 He atoms ad¬ 
sorbed on the upper side of ct-graphyne for upper-layer areal densities of (a) 0.0706 A -2 and (b) 
0.0941 A ~ 2 (red: high, blue: low). The black dots correspond to the carbon atoms and the white 
stars represent the peak positions of the lower-layer 4 He density distribution, which form a C 4/3 
commensurate solid. The computations were done at T = 0.5 K with the LJ substrate potential 
and the length unit is A. 
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FIG. 3: (Color online) Contour plot of two-dimensional density distribution of the upper 4 He 
layer at the areal density of 0.0687 A^ 2 , which corresponds to one less 4 He atoms per our 3x2 
rectangular simulation cell than the Mott-insulating density (red: high, blue: low). The black 
dots correspond to the carbon atoms and the white stars represent the density peaks of the lower- 
layer C 4/ / 3 commensurate solid. The computations were done at T = 0.5 K with the LJ substrate 
potential and the length unit is A. 



Z[A] 

FIG. 4: (Color online) One-dimensional density distributions of 4 He atoms adsorbed on both sides 
of a-graphyne as a function of the vertical coordinate z perpendicular to the graphyne surface, 
which were computed with the Yukawa-6 substrate potential. Here N up and N^ n represent the 
number of 4 He atoms per 3x2 simulation cell in the upper and the lower 4 He layer, respectively. 
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FIG. 5: (Color online) (a) Contour plots of two-dimensional density distributions of 4 He atoms 
adsorbed on the upper side of ct-graphyne for upper-layer areal densities of (a) 0.0706 A -2 and 
(b) 0.0941 A -2 (red: high, blue: low). The black dots correspond to the carbon atoms and the 
white stars represent the peak positions of the 4 He density distribution of the lower layer whose 
density is 0.0941 A" 2 for both cases. The computations were done at T = 0.5 K with the Yukawa-6 
substrate potential and the length unit is A. 
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FIG. 6 : (Color online) Contour plot of two-dimensional density distribution of the upper 4 He 
layer at the areal density of 0.0687 A^ 2 , which corresponds to one less 4 He atoms per our 3x2 
rectangular simulation cell than the Mott-insulating density (red: high, blue: low). The black dots 
correspond to the carbon atoms and the white stars represent the density peaks of the lower-layer 
C 4/3 commensurate solid. The computations were done at T = 0.5 K with the Yukawa -6 substrate 
potential and the length unit is A. 
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